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Abstract 

In this paper we construct a one-parametric family of (l+l)-dimensional one-component scalar field 
theory models supporting kinks. Inspired by the sine-Gordon and 4> 4 models, we look at all possible 
t— h extensions such that the kink second-order fluctuation operators are Schrodinger differential operators 

<H with Poschl- Teller potential wells. In this situation, the associated spectral problem is solvable and 

■^j- therefore we shall succeed in analyzing the kink stability completely and in computing the one-loop 

t-H quantum correction to the kink mass exactly. When the parameter is a natural number, the family 

becomes the hierarchy for which the potential wells are reflectionless, the two first levels of the hierarchy 
(— I being the sine-Gordon and 4 models. 

<D 1 Introduction 

From the very beginning of the surge of solitons in Quantum Field Theory in the mid seventies it 
was noticed that the second-order kink fluctuation operators in the sine-Gordon and ^ 4 models were 
Q\ respectively the two first members in the hierarchy of Schrodinger operators with "transparent" Poschl- 
Teller wells: 

O 

• K {1) = -— , + l-2sech 2 x , = -—. , + 4 - 6 sech 2 x 

■£j dx z dx l 

O 

In 1975, Christ and Lee in Reference [I] addressed the following natural problem: what should the 
potential energy densities in the parent field theoretical models be such that the kink second-order 
fluctuation operators are the other members in the "transparent" PT hierarchy: = — -pj + iV 2 — 

^ N(N + 1) sech 2 x, N G N? These authors found the general answer and described some analytical 
properties of the (/>^-kinks and the ((^-potentials arising in this kind of model. They also noticed 
that the parent potentials had branch points at the vacua for N > 3. In 1987, Trullinger and Flesch 
[2] succeeded in constructing the N £ N hierarchy by identifying the parent potentials in terms of some 
specific special functions. These functions are well defined when the field takes values in the interval in 
the field space between the two vacua. Beyond those points they tried to build the parent potentials by 
sewing the special functions just mentioned to some well-behaved functions with suitable properties at 
higher than the vacuum values of the field. The connection between the two zones, however, was at most 
twice differentiable at the branching points. Nevertheless, they showed that keeping only lowest-order 
(Gaussian) field fluctuations around the static vacua and kink solutions, the parent potentials are well 
behaved. This left open the question of computing the one-loop mass quantum corrections to the kinks 
associated with these parent potentials. In 1989 Boya and Casahorran applied the procedure of Cahill, 
Comtet and Glauber [3] to the models in this hierarchy. Thus, they computed the one-loop kink mass 
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quantum correction in terms of the bound state eigenvalues of the operator, see [4]. We stress that 

this formula is valid only if the potential well is reflectionless and suitably integrable. 

In a almost parallel development, the problem of the computation of the one-loop mass shifts for the 
sine-Gordon soliton and 4> kink was addressed. The first to tackle this were Dashen, Hasslacher and 
Neveu in 1974 in the framework of the ^-expansion of these (l+l)-dimensional field theories [5, 6, 7]. 
After a zero-point renormalization, which is performed by means of a mode-by-mode subtraction of 
the eigenvalues of the second-order fluctuation operators valued respectively on the kink and vacuum 
solutions and a subsequent mass renormalization, the authors derived the so-called first DHN formula. 
We stress that the DHN formula is again only valid when the kink second-order fluctuation operator is 
a reflectionless differential operator. In this case, both the kink and vacuum Hessians have a half-bound 
state whose contributions to the quantum correction annihilate each other in the mode-number cut-off 
regularization. In reference [8], the authors construct a generalized DHN formula which works smoothly 
when applied to quite general Schodinger operators. At the end of the past century a revival in the interest 
in this problem took place mainly due to role played by quantum solitons in supersymmetric theories. 
Many subtleties of the regularization/renormalization procedures involved were clarified, starting with 
the papers [9, 10, 11, 12]. Other related works, [13, 14, 15], addressed mixed issues in the problem by 
using different types of boundary conditions - PBC, Dirichlet, Robin-, regularization methods -energy 
cutoff, mode-number cutoff, high-derivatives-, and/or performed phase-shift analysis, in connection with 
possible modifications due to the quantum effects of the central charge of the SUSY algebra [16, 17]. 
In parallel, topological defect quantum fluctuations have been studied within the general framework of 
quantum vacuum fluctuations in recent years and the application of spectral heat kernel/zeta function 
regularization methods in this kind of problems has grown very effectively, see [18]. 

In the present paper, however, we address a two-fold goal: 

1. We extend the analysis to the case where the parameter a of the Poschl- Teller well is a positive 
real number, rather than a natural number, 

= — f^r + a 2 - a{a + 1) sech 2 x (1) 
ax z 

descending to the Christ-Lee- Trullinger-Flesch-Boya-Casahorran hierarchy when a = N £ N. We 
shall find analytically the kink solutions, which satisfy intriguing recurrence relations, in the form 
of regularized incomplete Beta functions. 

2. We compute the one-loop kink mass shifts by a variety of methods, depending on the difficulty of 
the problem. For a = N £ N the Cahill-Comtet-Glauber (CCG) formula [ ] is applicable but we 
shall also compute the quantum correction by the Dashen-Hasslacher-Neveu (DHN) mode-number 
cut-off regularized formula [5] and/or the exact heat kernel/zeta function regularization method 
[19]. If a is not an integer, only the generalized DHN procedure as explained in [8] can be applied. 
Alternatively, the asymptotic expansion of the heat kernel will allow us to obtain sufficiently good 
approximations to the one-loop kink mass quantum correction. 

The organization of the paper is as follows: in Section §. 2 we shall construct the family of models 
and we shall also identify the different kinks, discussing interesting relationships between the kinks in 
different models of the family. We shall also study kink stability, which by construction reduces to the 
spectral analysis of the operator (1). In Section §. 3 we shall compute the one-loop quantum correction 
to the kink masses. We shall distinguish the cases a £ N and a ^ N, or equivalently, we shall deal with 
the Hessian operator (1) with reflection and reflectionless potential well separately. We shall approach 
this calculation from different points of view. In particular, we shall use the CCG (Cahill, Comtet 
and Glauber) formula, the DHN (Dashen, Hasslacher and Neveu) formula, the spectral zeta function 
regularization method and the asymptotic approach. The first three methods can only be applied to 
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models that involve reflectionless Hessian operators, although the second one has been generalized in 
reference [8]. A comparison between the results obtained with the different methods will be offered, 
showing a high level of accuracy. In the last Section we shall summarize our conclusions. 



2 The family of models: kinks and their stability 

2.1 Generalities 

The action governing the dynamics in our (l+l)-dimensional relativistic one-scalar field theoretical mod- 
els is of the form . 

sw = //*%■ 

Here, ipiu 11 ) '■ ^ -> R is a real scalar field; i.e., a continuous map from the (1 + l)-dimensional 
Minkowski space-time to the field of the real numbers. y° = r and y 1 = y are local coordinates in 
IR 1 ' 1 , which is equipped with a metric tensor g^ u = diag(l, — 1), fj,, u = 0,1 such that y^ = g^y^yu, 

d d fj,i> d d 

We shall work in a system of units where the speed of light is set to one, c = 1, but we shall keep 
the Planck constant h explicit because we shall search for one-loop corrections, proportional to H, to the 
classical kink masses. In this system, the physical dimensions are: 

[h] = [S]=ML , [ Vli ]=L , [ijj]=MH 1 2 , [U] = ML- 1 . 

The models that we shall consider are distinguished by different choices of the part of the potential 
energy density which is independent on the field spatial derivatives: U[ip(y^)]. In all of them there will 
be two special parameters, and 7,2, to be determined in each case, carrying the physical dimensions: 
[rrid] = L^ 1 and [7J = M~ %L~ 2. We define the non-dimensional coordinates, fields and potential in 
terms of these parameters: 

Xn = m d y ti , x = t , xt = x , <j) = -y d tp , U(<p) = -%U(tf) 

m d 

The action S[ip\ is also proportional to a non-dimensional action, i.e., S[ip] = \S[cf>], where 

Id 



S[<j)} = / dx 



l -d^d^-U{^) 



(2) 



Here 4>{xp) : R 1 ' — > K is a non-dimensional real scalar. Also, we define a non-dimensional potential 
energy functional by means of the relationship E[<j>], where 



E[4>] = 3*-EW\ , E[4>] = r dxe(x) = r dx \ (®f \ + U{4>) 



(3) 



From now on we shall work in this non-dimensional context. We shall assume that U(4>) is a non-negative 
twice continuously differentiable function of <j) G K, i.e., U{(j)) G C 2 (IR) and U((f>) > 0. We now want to 
construct kink supporting models where the kink stability is governed by Poschl- Teller spectral problems. 
We shall thus assume that the U(4>) has at least two minima -which can be fixed as the zeroes of U((j))- 
at the values <f> = ±1. This assumption does not imply a loss of generality because a translation of the 
origin and a dilatation in the field space allows us to set this prescription in an arbitrary model. 

The conventions are those set out in reference [20]. 
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The configuration space C of the system is the set of the maps from the spatial line K to the in- 
ternal (field) space R for a fixed time t, such that the energy functional E[<f>] is finite, C = {4>(t, x) S 
Maps(M 1 ,IR) : E[(f>] < oo}. The search for kinks is tantamount to the search for finite energy static 
solutions <p(x) of the field equations, which in this case become the "Newton" ODE: 

= dU 

dx 2 d<f> { ' 

A Lorentz transformation shows the static solution 4>(x) in its center of mass as the traveling wave: 
<p(t,x) = 4>{ )■ Moreover, spatial translations and reflections, x — > x + xq, x \-¥ —x, provide more 
solutions because of the invariance of (4) with respect to these transformations. Therefore, the complete 
family of solitary wave solutions of (4) is compactly written as: 4>(x), x = (—l) a x ^°~% t with a = 0, 1. 



The second-order ODE (4) always admits the "mechanical energy" I\ = ^ — U{<f>) as a first 

integral. Because the "mechanical action" is the field theoretical energy, the energy finiteness requirement 
identifies the kink solutions as the trajectories of I\ = mechanical energy: 

f x = ±vmf) ■ (5) 

Therefore, the kink solutions interpolate between the (fry = =Fl vacuum at X = — oo and the <j)y = ±1 
vacuum at x = +oo: 

lim (pxit, x) = =Fl e M , lim (j) K (t,x) = ±1 e M 

From the signs in the equation (5) one sees that the interpolation is monotonic: the kink solutions 
are either monotonically increasing (kinks) or decreasing (anti-kinks) functions of x from =pl at the left 
boundary of the line to ±1 on the opposite spatial limit. 

The linear stability analysis of any static solution - vacua (fry or kink ipgix)- is codified by the spec- 
tral problems Kof UJO (x) = uj^f^^x) and Kf LJ (x) = co 2 fu{x) for the associated second-order fluctuation 
(Hessian) operators: 

d 2 d 2 U ri . r , d 2 d 2 U ri , N1 
K ° = -dx- 2 + W [M ' K = -dx- 2 + W^ K[x)] • (6) 

The stability of the solutions 4>y and 4>k(x) is guaranteed when the eigenvalues of these operators are 
non-negative. In particular, we are interested in the "parent" potentials that admit kink solutions such 
that: 

= -^+a 2 - <t(<7 + 1) sech 2 x , (7) 
dx z 

where a G M. The potential wells associated with these operators V^ix) = —a{a+\) sech 2 x are depicted 
in Figure 1 where the integer values of a has been emphasized. Note that the wells become deeper and 
higher as a increases. 

We shall denote by c/)^Z\x) the kink solutions and by U( a \(j)) the parent potential for a given a in 

this family of models. Because the kink solutions asymptotically reach the vacuum values 4>y^ = ±1, the 
kink Hessian operators are accompanied by the vacuum Hessians: 

K (a) --— +a 2 
K ° ~ dx 2+ ' 

Thus, the mass of the fundamental quanta in these models is: %M-[4>v] = ■ 
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Figure 1: Graphical representations of: (a) the potential well in the Hessian operator for several values of a (b) 
the kink profiles for several values of a. The cases a — N have been emphasized. 



2.2 Zero modes and kinks 

Our first task is the identification of the parent potentials such that: 

[<j$ (x)] = a 2 - a{a + 1) sech 2 x . (8) 

The solution of (8) seems hopeless without knowledge of <jr^ (x). The equation (5) 

^ = , (9) 

however, does allow us to interpret (8) as an ODE in the x variable via the change of variables from 

the implicit parent potential for the kink solutions by the explicit function of x il a (x) = (cj)^\x)) 

7W^~dr d<f> 2 ~ Wv^d^ ~ (2U CT ) 2 



to x. Recall that <p"p : IR — > [—1, 1], x \— > (f>^\x) is a bijective map determined by equation (9). Trading 
the implicit parent potential for the kink solutions by the explicit function of x il a (x) = (cj)^\x)) 
and bearing in mind from equation (5) that ^ = and ^ = ^<gf - j^i^) 2 , the 

equation (8) becomes: 

i d 2 ii a : fdii a \ 2 a 2_ a{a + 1)sech 2- (10) 



2iX a {x) dx 2 (2H (T (x)) 2 V dx 



(8) is a non-linear ordinary differential equation that can easily be reduced to the Riccati equation by 
means of the change F{x) = \ logil CT : 

d 2 F (dF\ 2 2 , l2 _ 
^ + =a -a(a + l)sechx . 

A particular solution can be obtained by the standard procedure of sending the Riccati equation to the 
Schrddinger equation via the new change 4^ = dlc f x ^ , or, il = -ip 2 .: 

= i ' 2 ~ a ( cr + l)sech 2 x) ipo , (11) 

i.e., ipo(x) = \J C a sech. a (x) is the ground state of zero energy in the PT potential well. Therefore, we find 
a solution depending on the integration constant C a of our problem: 

tt a (x) = Casech^x , Jim U tT (x) = Q . (12) 

x— >±oo 

This solution fits smoothly into our scheme because the kinks tend to the zeroes of the potential (</>) 
at x = ±oo. The hidden reason for finding this expression (12) for the parent potential as a function of 
the parameter x is that via the equation 

#£ } = y/2i^{x)dx (13) 
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the potential is related to the translational mode J£ , which is in turn the ground state of zero energy 
of KM because of general arguments (just plug it in (4)). Integration of (13) immediately produces the 
kink solution 

= J^ T [|] tanh x 2 Fi[|,l-f,f, tanh 2 x] (14) 
in terms of a Gauss hypergeometric function. The constant C a has been tuned to 



2r[| + § ] 2 



7rr[f]2 

in order to set the zeroes of the parent potential at the <j) = ±1 values. One can check that for a = 1 the 
soliton of the sine-Gordon model (j)^(x) = ^ arctantanh | is found, whereas the kink of the </> 4 model 

(2) 

(jr K (x) = tanhx arises when a = 2. 

It is convenient to rewrite the kink profile (14) using other special functions. We recall the definition of 
the beta function in terms of three Gauss Gamma functions: B(a, b) = t a ~ l (l — t) b ~ l dt = r fa+ffi • We- 
also need the incomplete beta function B(z; a, b) of a free upper integration limit z € [0, 1]: B(z; a, b) = 
j a: t a ~ 1 (l — t) b ~ 1 dt = \z a 2-Fi[a, 1 — 6, a + 1, z], which is related to the hypergeometric Gauss function. 

Finally, we shall consider the regularized incomplete beta function I(z; a, b) = ~g7^gr ■ With this arsenal 
we obtain the very compact formula: 

(x) = sign(x) /[tanh 2 x; \ , f ] . (15) 

Note that formula (15) is an odd function because of the sign function that multiplies I[tanh 2 x; |, f]- 
This does not induce singularities because / is zero at the origin, see Figure 2. It is remarkable that the 
soliton and kink solutions are embraced in this general formula. 

2.3 Recurrence relations 



The magic of the hypergeometric functions gives rise to interesting recurrence relations between kinks for 

z d 
b dz ' 



different values of a from (14). The well-known formula 2-F1 [a, 6+1, c, z\ = 2-F1 [a, b, c, z] + f ^2^1 [a, b, c, z] 
applied in (14) means that: 



(f> ( ° +2) (x) = #(z) + - tanh a; ^-(x) , (16) 
a dx 

a relation between the a and a + 2 kink profiles. Taking into account (13) we rewrite (16) as: 

^ +2 \x) = 4%) + 2r fe| tanhx sech-x , (17) 

a form of the recurrence relations that can also be derived from the identities between regularized beta 
functions: I(z,a,b) = 1 — 1(1 — z,b,a) and I(z,a + 1,6) = I(z,a,b) — , applied in the formula 

(15). In formula (17) one notices that the (a + 2)-kink is equal to the addition of the a- kink plus an odd 
function, implying that | </>^ +2 ^ (x) \ > \4>^\x)\. Therefore, the kink profile is more and more concentrated 
around its center in members of the family with higher and higher a, see Figure 1(b). 

For odd and even integer values of a the formulas simplify analytically. We discuss the a = 2k even 
and a = 2k + 1 odd, k £ N, cases separately. If a = 2k, (17) becomes: 

0f + 2 ) (x) = 4 k) (x) + {2k 2 ; k f A tanh x sech 2fc x . 
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In the odd case the same recurrence relation reads: 



{2fc+1) (x) + ^±^tanhxsech 2 ^3 



In Table 1 we show the kink profiles of the first ten members of the hierarchy of models. 



%\x) 



- arcsin tanh x 

7T 

0# 0*0 + ~ tanh x sech x 
<pj^{x) + jp tanhxsech 3 x 
0^(x) + tanhxsech 5 x 
<Pk\x) + tanh x sech 7 x 





<4%) 


2 




tanhx 




4 


4\x) = 


4\x) 


+ 2 tanh x sech 2 x 


6 






+ | tanh x sech 4 x 


8 


4\*) = 




+ tanh x sech 6 x 


10 


ff\x) - 


-4\x 


) + jf^ tanhxsech 8 x 



Table 1: Some kink solutions of the hierarchy of models characterized by a G N. 

There also exist recurrence relations between a- and a + 1-kinks differing in one unit in the family 
parameter, although in this case the formula involves a quadrature. From (13) we obtain 



\AWi(a 



a 



r[ 



2 m + 



£12 

21 



£12 



sech x i 



which can be two-sided integrated to find: 



a 



r[f] 2 
2r[i + fp 



sech x </>^ (x) + / sech x tanh x (p^ , (x)dx 



An integration by parts has been performed in the right-hand side member. 

The classical energy of the kinks can easily be worked out from the zero modes, the connection 
formulas between hypergeometric functions and the Kummer theorem: 







LI M ( 




dx ) 





4 <r+i r[| + §] 2 



• 2 



Fi[(T, 2<t,1 + ct;-1] 



4r[q]r[i + f ] 2 
0Fr[a + i]r[f] 2 



2.4 Parent potentials 

To completely define the family of models we identify the potential (12) as a function of the field 
first, inverting formula (15): 



x 



1.1 £ 

J K I' 2' 2 



arctanh W / 1 

and, then, plugging this expression into the formula (12) to find: 



vrrff] 2 



by, 



(18) 



We emphasize that this formula is only valid in the interval <j) G [—1,1] in the internal space, see Figure 
2(a). The reason is that the kink solutions live only in this interval and use of the map 4$ : R — > [—1, 1] 
to reconstruct the potential in (2) is restricted to this interval. Obviously the character of the potential 
for field configurations such that \4>\ > 1 does not influence the kink solutions. We have the freedom, 
however, to mathematically extend the definition of the potential as a function of to the \(f)\ > 1 
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range. Even requiring twice-differentiable functions U^ a \(f)) E C 2 (M), there are several possible choices 
of continuations of U^ a \(j)): see [2] for some examples. In Figure 2(b) we show periodic continuations of 
the potential defined in <j) € [—1, 1] to the \(j)\ > 1 half-lines. Here, the piece of U depicted in the Figure 
2(a) is cloned in all the [2k — 1, 2k + 1], k E Z, intervals. The continuations depicted in Figure 2(c) are 

accomplished analytically by expanding '((/)) = ^[^f 2 ~~ 0^1' 2' f as a P ower series near 
the (/> = ±1 points and truncating the series to a polynomial at a finite higher than two order. 




Figure 2: Graphical representation for the potential: (a) in the range \(f>\ < 1, (b) periodic continuations, (c) 
polynomial continuations, (d) a = 4 generalized Chebychev continuation, to \<j)\ > 1. 



Starting from the a = 1 potential in formula (18) extended by the periodic continuation scheme we 
obtain: U^'((j)) = + cos(n<f>)], a beautiful version of the sine-Gordon model (properly translated 
in internal space). Extending the a = 2 potential in (18) by using the second type of continuation we 
recover the (/> 4 model: U^ 2 \(j)) = \{4> 2 — l) 2 . Indeed, these two examples are the exceptional C°°(M) 
continuations in the whole hierarchy of a = N potentials. More generally, C 3 (M) continuations of the 
potential (18) are not possible if a > 2 because 



d 3 U^ 



d 



[a 2 — a (a + 1) sech 2 x] = 2a(a + 1) sech 2 x tanh x 



_dx 2a(a + 1) 



sech 2 ""x tanh x 



diverges when x tends to infinity (the vacuum points 4> = ±1). There is however another case, a = 4, 



where the (18) U^(4>) = § [l - 2cos(| arcsin < 
in terms of elementary functions: 



extended to the whole real line K can be expressed 



1 — 2 cos (| arcsin 0)] 4 if \4>\ < 1 
1 — 2 cos ( | arcsin ^ J if \<f>\ > 1 



(19) 



This function, akin to Chebychev polynomials, is depicted in Figure 2(d) and is easily seen to be C 2 (R). 
Note, however, that this continuation differs from that depicted in Figure 2(c), in particular the extended 
potential tends to f at <j> = ±oo. 
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2.5 Kink stability 



We recall that the Hessian differential operator at the vacuum solutions <py ' = ±1 reads 



K, 



d 2 c>W 



dx 2 



+ 



[±1] 



dx 2 



see (6). In order to analyze the spectrum of this differential operator we temporarily confine the system 
in an interval of (non dimensional) length I and impose periodic boundary conditions on the plane eigen- 



waves i(j^(x) = ~^\ e %kX • The wave vectors are quantized kl = 2irn,n 6 Z, and in the limit of very 

large I the spectral density is: p(K^ ) = = The corresponding eigenvalues are co 2 (k) = k 2 + a 2 
such that the threshold of the continuous spectrum arising at I = oo is a 2 . The constant function is 
also a non-propagating eigen- function of Kq with energy a 2 . Bound states buried in the threshold of 
the continuous spectrum are called half-bound states because due to the ID Levinson Theorem must be 
counted with a weight of one-half. In sum, we have that Spec(K { Q a) ) = {a 2 }i U {k 2 + a 2 } km . All the 
eigenvalues are positive, and there are no tachyonic waves fluctuating around the vacua, which accordingly 
are stable. 

The analysis of kink stability requires us to deal with a more complicated but tractable spectral 
problem: K^ip LJ (x) = uP^^ix). We confront the Schrodinger operator 



K^ = ~ + v* + V{x) 



dx 2 



+ a 2 — a (a + 1) sech 2 £ 



(20) 



li 2 



i.e., the quantum Hamiltonian describing the motion of a particle of (non-dimensional) mass fj, = in a 
x — > —x symmetric Poschl- Teller well (7). The scattering eigen- functions of eigenvalue ui 2 {k;a) = k 2 + a 2 
are well known [21]: 



^\x) = 2 ik C { k a) sech- ik x ■ 2 F X 



a + 1 — ik, —a — ik,l — ik, 



k g 



(21) 



These eigen-functions belonging to the continuous spectrum of can also be written in the simpler 
form 



^\x) = ci a) e ikx 2 F 1 



-a, a + 1, 1 — ik) 



e x -\-e~ x 



(22) 



by using the linear transformation: %F\ [a, b, c, z] = (1 — z) c a b 2F\ [c — a,c — b, c, z] . From the asymptotic 
behaviour 



~ (~i{<t) ikx 

X— ¥00 ^f. c 



-x— >— OO 



r(l - ik)T(-ik) 



e ikx , T(l ik)T(ik) ^jkx c -ikx 



T(a + 1 - ik)T(-a - ik) r(cr + l)r(-<7 

read in (22), the transmission and reflection amplitudes are identified: 

t (*)r k \ _ r(<r + 1 - zfc)r(-<T - ifc) r («r)(jfe) = r (°" + 1 - ik)T(-a - ik)T(ik) 



r(l - ik)T{-ik) 



r(l + a)T{-a)Y{-ik) 



(23) 



The phase shifts in the even, 5+'(k), and odd, S_(k), channels are related to these amplitudes via the 
identities: 



e 2i5 ±\k) = t^±A a \k), 5 { ±\k) = ^arctan 



1m(tW ±A a \k)) 
Re^M ±r( ff )(fc)) 

t( ff )(jfe) = cos (S^{k) - 5 < f \k)j ■ exp[UsP(k) + <5 M (fc))] 
A a \k) = sm^\k)-5 W) (k)yexp[^(5^\k)+5 ( L\k)+Tr)] . 
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Again we impose PBC e ikx = e ik{x+l ^ ■ e £(*+ } (*)+*- 3 (*)) such that kl + \{5 { + ] \k) + 5 W (A;)) = 2vm and 
the spectral density of the kink fluctuations is, for very large I: p(K^) = 4jr = ^- + ^^jf-(k), where 
5^ a \k) = 5+\k) + 5_(k) are the total phase shifts -the sums of the phase shifts in the even and odd 
channels- of the scattering eigenfunctions of . 

Strictly at the I = oo limit the spectrum of includes a discrete part. There are several bound 
states: k = j{a — j),j = 0, 1, 2, • • • , I [a], are I [a] + 1 poles of t^ a \k) in the positive imaginary axis of 
the complex A;-plane. Here, I [a] denotes the integer part of a and we temporarily assume that I [a] ^ a. 
uij(a) = —(a — j) 2 + a 2 = j(2a — j) are the eigenvalues of these I [a] + 1 bound states, whereas the 
eigenfunctions, derived from (21) for a better visualization of the L 2 -integrability character, read: 

tf\x) = {eX 2F,[2a + 1 - j, -j, a + . 

Note that u^(ct) = 0, Vex. The ever present zero-mode eigen-function is the Goldstone boson owing to the 
spontaneous breaking of the translational symmetry by the kink. 

In the special cases when a = N G N is a natural number two related novelties arise: 

1. The reflection amplitude r N \k) = becomes null and the kink potential well is transparent to the 
motion of mesons. The transition amplitude, however, becomes: 

fWm n iv-i N-j-ik _ N _! k 2 -(N- j) 2 + 2ik(N - j) 

t W- il i=0 N _ j + ik - U -j=Q {N _ j) 2 + k 2 ■ W 

2. There are only bound states in the spectrum of K^ N \ although a half-bound state just buried 
at the threshold A 2 of the continuous spectrum arises: 



^ JVj (5:)=Cr2^l[A + l,-iV,l, 

2 2 



In sum, the spectrum of is slightly different in two situations: (1) if a = N £ N we have: 
Spec(AW) = {j(2N - j)} j=0 ,i,...,N-i U {A 2 }i U {k 2 + N 2 } km . (2) if a $ N, Spec(K^) = {j(2a - 

j)}j=o,l,.„,l[a-] U {k 2 + °" 2 }A;eiR- Therefore, the spectrum of the Hessian operator is non-negative and 
the kink solutions are stable in all these models. 



3 One-loop quantum correction to the kink masses 

In this Section we address the main issue in this work: computation of the quantum correction to the 
classical kink masses for arbitrary values of a G M + . We stress that as far as one- loop (Gaussians) 
fluctuations are concerned this is a well defined problem, see [1, 2, 4]. It seems to us convenient to 
distinguish the a = N G N and a ^ N cases in these computations. For the sake of completeness and 
to gain better knowledge of this subtle quantum phenomenon we shall address this task using different 
tools: the CCG formula, the DHN formula, the generalized zeta function regularization method, and the 
asymptotic heat kernel approach. The four methods lead to compatible results, reinforcing the reliability 
of each individual method and the conclusions obtained. 

The one-loop kink mass -the classical kink mass plus the one-loop quantum correction- is of the form: 

£[V^]i-ioo P = ^ {®\4>k\ + A • 
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Thus, formally, the object to be computed is: 



AM 



:Tr 



L 2 



(K 



K 2 



where the operators entering formula (25) are: 



-Tr 



L 2 



(25) 



dx 2 dd> 2 



dx 2 



K 



dx 2 



+ 



d 2 u 



4>K 



dx 2 



+ v 2 + V{x) . (26) 



We must make sense of the L 2 -trace of quotients and differences of square roots of differential operators 
to account for the effect of the kink quantum fluctuations measured with respect to the vacuum quantum 
fluctuations. In order to do so it is convenient to divide AE[<fyjg] = AEi[(j)^ ] + AE^^x ] ^ wo P ar ts 
with different physical origins: 



1. AE 1 [4 ) 



^Tt L 2(K2 — Kq) is due to the zero-point renormalization of the kink vacuum energy 
by subtracting the vacuum energy. We emphasize that this formula means that the subtraction is 



performed mode- by- mode. In infinite dimensional spaces Ti L 2(K2 — Kq) ^ Tr L 2Kz — Tr L 2K 2 ; 



the difference in the traces (not the trace of the difference) would arise if the total energy of the 
vacuum fluctuations were subtracted from the total energy of the kink vacuum fluctuations. This 
is a source of many controversies in the Literature. 



2. A£ 2 [#] 



lTr L2 (^¥ 



comes from the one-loop counter-term induced by the mass renormal- 



ization procedure in (1 + 1)D scalar field theory merely due to normal-ordering of the Hamiltonian. 



3.1 The a = N G N case 

3.1.1 The Cahill-Comtet-Glauber (CCG) formula 

We shall always deal with potentials such that: dx (1 + |x|)y(x) < +oo. If a = N £ N, the reflection 
amplitude is null, A N \k) = and the potential well is transparent. Given these two circumstances, it is 
shown in Reference [3], relying on Zakharov and Faddeev, [22], that: 

1. There is a nice relation between the bound state eigenvalues and the phase shifts given by the 
formula: 

6j = arccos^- . (27) 
v 

j 

Here, uj 2 are the bound state eigenvalues of K, i.e., those below the threshold of the continuous 
spectrum: < oj 2 < v 2 . The sum is, of course, over all these eigenvalues. 

2. The total area enclosed by V(x) is also related to the bound state eigenvalues: 

/oo 
dx V(x) = -4v sin Oj . (28) 
-oo ■ 



j(*)(Jfc) =2^ arctan 



.k 



sin ( 



3. There is a last feature in the spectrum of these reflectionless potentials to be taken into account: the 
highest bound state of K sits exactly at the threshold of the continuous spectrum. The contribution 
of these eigen-functions -half-bound states- cancel exactly with the contribution of the constant 
eigenfunctions -also half-bound states- of Kq. 
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With all this information at our disposal, we derive the CCG formula. The first contribution AE\[(j) K ' 
is 



1 il -I roo 

-Tr„(K*-K§) = - £ Uj - -J ^ dkfr\k) 
— lim < ( \/k 2 + u 2 arctan[— sin Q~ 

>7T A^oo^ [V l k J 



V 



E 



COS I 



2tt 

7T 



E[ 



arctan[cot#j] J cos9j — sinOj 



v 

2^ 



J dk a/A; 2 + v 2 — ^arctanf- sin 1 

Esin^i lim log 
J A->oo 



A + VA 2 + 
-A + VA 2 + t; 2 



while the second AE^lfi^] reads 



V f V ( x ) 
-Tr i2 ^ — 



K, 



1 



dx V{x) 



dk 



o 



2tt 



Esin^i lim log 
J A^oo 



A + VA 2 + v 2 
-A + VA 2 + v 2 



E 



sin ( 



dk 



Vk 2 + v 2 



Adding these two pieces together the logarithmic divergences cancel and we obtain the very compact 
formula 

= > (sin i 

7T " 

3 



3 

for the one-loop kink mass shift derived in [3] in terms of the threshold v and the bound state eigenvalues 



4 



In the models studied in this paper the operator is of transparent Poschl- Teller type, supporting 
N bound states with CCG angles: 



9j = arccos 



Vj(2N-j) 
N 



Recall that u](N) = j(2N-j) and v 2 = N 2 . Therefore, formula (29), after some manipulations, provides 
the exact result in the form: 



AM 



N 



TT 



N-l 

E 

j=0 



sin arccos 



N 



N-l 



-—N(N + 1) + - V" \/ N 2 — r 2 arcsin 



r=l 



N 
r 

N 



arccos 



N 



(30) 



where we have defined r = N — j . 

The very well known results for the sine-Gordon- and </> 4 kinks are reproduced here. If N = 1, we 
find the one-loop quantum correction to the sG-soliton mass: AE[<p^] = — i. For N = 2, we obtain the 
one-loop quantum correction to the c/> 4 -kink mass: AE\<f>g\ = — ^ (in this case = but in 
general rrid = 111 Table 2 we display the results obtained by means of (30) from N = 1 to N = 20. 

A graphical representation of the dependence of the one-loop mass shift in N is also depicted in the next 
Figure. 
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4 


-1.58003 




14 


-11.0575 


5 


-2.15555 




15 


-12.4651 


6 


-2.81245 




16 


-13.9568 


7 


-3.55156 




17 


-15.5328 


8 


-4.37342 




18 


-17.1929 


9 


-5.27839 




19 


-18.9374 


10 


-6.26675 




20 


-20.7661 



Table 2: Numerical values of the one-loop kink mass quantum corrections (15), AE[<f>^] for N = 1, . . 
graphical representation of the one-loop kink mass shift as a function of N. 



, 20 and 



3.1.2 The DHN formula 

In the seminal series of papers [5, 6, 7] Dashen, Hasslacher and Neveu computed the one-loop quantum 
correction to the kink and soliton masses by using what is known as the "first" DHN formula 2 . Recently, 
the DHN prescription has been thoroughly revised by several research groups in order to address similar 
issues about supersymmetric kinks [9, 11, 13, 14, 17]. In this framework several subtleties arise, mainly, 
but not only, concerned with compatibilities between boundary conditions, regularization procedures, 
and supersymmetry. Like the CCG formula, which is no more than an elegant reformulation of the 
DHN result in terms of the bound state eigenvalues and the threshold of the continuous spectrum, the 
DHN procedure works fine when the potential wells involved are transparent. In this case, the DHN 
expresses the same L 2 -trace combinations without using the connection between phase shifts and bound 
state eigenvalues encoded in the CCG formula: 

AE[<f> K ] = lTr i2 (id-ir|)-Tr i2 ^4= (31) 

Because the kink Poschl- Teller wells are transparent when a = N, the phase shifts 5^ N \k) are only 
identified from the transmission amplitudes t^ N \k) (24) and, in turn, the kink spectral density is easily 
derived: 



1 JW 2vr ^ 2vr dk 2ir ir ^ r 2 + k 2 

r—l 



The area of the kink well (x) = -N(N+ 1) sech 2 x is also easy to compute: (V^ N \x)) = -2N(N+1). 
Use of the identity N (^ +1 ) = YlrLi r anc ^ cancelation of the r = N terms allow the calculation of the 
integrals involved in formula (31): 

dk- — — \/ k + N 2 — — (V(x)) 



2ir J 2vr dk 8tt ' y " ^k 2 + N 2 




N 

y r 

^— ' k 2 + r 2 

r=l 



2Vt 2 + JV 2 



The "second" DHN formula offered WKB computations on the energies of quantum sine-Gordon breathers. 
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A N-1 
lim / dk > 

x ^ r (r 2 -iV 2 ) 



E 

r=l 



r 

7T 

lim 



VA: 2 + N 2 



+ 



1 



7T 



A— >c 



k 2 + r 2 Vfc 2 + iV 2 



(A 2 + r 2 )\Jk 2 + N 2 



N-1 

E 

r=l 



ViV 2 - r 2 



arccos 



7T 



AT 



Therefore, after adding the bound state contribution and subtracting the boundary term arising in the 
mode-by mode subtraction, in the DHN approach the one-loop quantum correction to the kink mass 
reads: 



N-i 



1 vOV 2 - 



i=l 

■ijV(iV + l) + ^^ 



2vr 

AT-l 



r=l 



arccos — 

7T TV 



r=l 
N-1 



tt r 

arccos — 

2 N 



—N(N + + y/jv 2 - r 2 arcsin 

27T 7T ^-f 



r=l 



r 
iV 



We have now performed the inverse change: j = N — r to find an identical result to that offered by the 
CCG formula (30). 



3.1.3 The spectral zeta function procedure 

Traces and determinants of powers of elliptic operators are only defined by means of a process of analytic 
continuation that mimics the definition of the Riemann zeta function as a meromorphic function, giving 
formal meaning to strictly divergent series by extending the usually rational exponents of the eigenvalues 
to the complex plane. Replacing natural numbers by the eigenvalues of positive elliptic differential 

— 2 

operators, the spectral or generalized zeta functions are defined. Assuming that A G A\L (R)] is an 
operator with a positive discrete spectrum, Spec A = {lo 2 £ R : Aip n = uj 2 ip n }, the spectral ^4-zeta 
function is the formal series obtained by summing over the — s-power of the eigenvalues: 

C A (s) = TvA~ s = Y,-^ , . (32) 

n ° Jn 

There is a general theory of elliptic pseudo-differential operators that characterizes the conditions under 
which the generalized zeta function is a meromorphic function. By a process of analytic continuation, 
close in spirit to the analytic continuation defining the Riemann zeta function, the values that are not 
at the poles of the spectral zeta functions and their derivatives are taken as "regularized" definitions of 
traces and logarithms of determinants of (complex powers of pseudo-)differential operators [23, 24, 25]. 

We closely follow this idea to regularize the kink Casimir energy by using the spectral zeta functions 
of the K and Kq operators: 

i — 00 \ 1 

AExfaM = -Tr L2 (V) S - Tr L2 AV = - £ W (n)- 2s - £ „ (n)- 2s = -[C K± (a) - (k (s)] 

\n=l n=l / 

(33) 

appart from the physical pole: s = — |, [19]. By K 1 - we denote the differential operator K acting in 
the sub-space of the Hilbert space orthogonal to its kernel. From a mathematical point of view, zero 
modes are disastrous in the definition of spectral zeta functions, giving rise to uncontrollable divergences. 
Fortunately, physical reasons show that zero modes due to the symmetry breaking of some continuous 
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group only enter in two- and higher- orders in the loop expansion, see e.g. [26]. Simili modo, the other 
contribution coming from normal-ordered Hamiltonians can be also regularized: 

AE 2 ^ K )[s] = -\Tr L2 [Vix)^- 1 ] = \ (V(x)) lira |E!j^ CjCo ( a + 1) . 

There is a very subtle issue here. By choosing the regularization point in s + 1 instead of s we pick 
the physical pole of both divergent quantities at the same point: s = — \ in C. A miracle happens: 
the different residua provide, after cancelation of the divergent parts, exactly the same correction as the 
mode-by-mode subtraction, not present in the (33) term. In sum, the one-loop kink mass shift is obtained 
by this procedure through the formula: 



AE(cf> K )= lim [A£!(fe)[ S ] + A£ 2 (fe)[ S ]] 



It is also convenient to use another spectral function, the i^-heat or partition function, because as a 
series in decaying exponentials of the eigenvalues this function enjoys better convergence properties and 
is related to the spectral zeta function via Mellin's transform: 



M/?) = Tr i2 e-^ = 5>-^(") , Ck(s) = -^ 



dpp s - l h K {p) , 



where f3 is a fictitious inverse temperature. In our models the JC-heat function can be identified analyti- 
cally from the spectral data of the kink well: 



Tr L2 e 



N 

E 

j=0 



e -/3j(2N-j) + 



I 

2^ 



N 



N 

E 



/ r 2 + k 2 

r=l 



-/3N 2 



+ ^ E rf [ry/0\ 



r=l 



as a sum of error functions Erf [z\. Because of the zero point vacuum renormalization and the exclusion 
of the zero mode the important quantity in the kink mass shift problem is: 



N 

h KW ,W)-h K(N) (P) = ^L^ KW -Tr L2 e-^ N) =e-^^e^Erf(rv^)-l 



(34) 



r=l 



a magnitude that is plotted in Figure 3 as a function of (3. The Mellin transform of formula (34) reads: 






Figure 3: Graphical representation of the difference between the K^ N ^> -heat and K^-Yie&t functions: (a) the 
exact function for N = 1,2,3,4,5, (b) the exact function for N — 4, (c) the asymptotic approximation for several 
N t truncation orders. Increasing values of N t give better approximations to the exact function. 
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JV— 1 



c K w±(s) - ck (s) = e 

r=l 
N-l 

= £ 



(N 2 — r 2 )~ s 



1 r[ S +i] , ^2 

2 Fi[s,s+ |,s + 1, 1 



0Fr 2s r[s + 1] 



ttN 2s sT\ 



r=l 



2s I _/y"2 _ r 2 



2+s r[ s + 



r[«] 



2J n fl , 1 3 H 1 

2-Tl[ 2) * i- 2 ) 2' AT2_ r 2j 



rf s + 



Bearing in mind that the spectral Ko-zet& function is £ (jv)(s 
one-loop shift due to mass renormalization becomes: 



2v^ 



-L- (N 2 )l s ^-fepp-, the regularized 



iV(iV+1) (iV 2 )-|- sr[S+ 2 ] 



2y^ V ' T[s] 
The Laurent expansions around the pole s = -\ of both AEx(<$ } )[s] and A^ 2 (^ ) )[s 

r 1 



Af-l 



A^l(C>] = E 



r=l 

iV 1 



2vr s + A 2tt 



T + 7T h + 1 °S(N 2 - r 2 ) + ^(-i) - 2 F[[\, 0, |, ' " 



7V 2 -r 2 ■ 



^ + r + — (-2 + 7 + 2 log iV + ^(-i)) + 0( a + i) 



N(N + 1) 1 
4vr s + i 



+ 



Ar(W+1) (7 + «-i))-^ + iiVlogiV= + 

7T Z7T 



47T 



Af-l 



Af-l 



+ ^ E no g (N 2 - r 2 ) - I E r • 2 i^,0, §, 



r=l 



r=l 



A^ 2 (4 Ar) )H 



pave the way to follow a cautious path falling in the physical limit: 



lim 

N-l 



AE 1 {^ ) )[s]+AE 2 {<j>^>)[s] 

N-l 

2tt 



r=l 



N 1 

7T 27T 



Nlo R N 2 - N{N + 1) lo,N 2 + 



4tt 



+^E H °g^ 2 -r 2 )-^E-^[io, 



2 ' N 2 -r 2 t 



(35) 



r=l 



Here, 2 Fl[a,b,c, z] stands for the derivative of the hypergeometric function with respect to the second 
argument; 7 is the Euler number, and ^{z) denotes the digamma function. Note that AE\{4>^)[s\ 

and A£' 2 ((/)^ V )[s] have a simple pole at s = — \ with opposite strengths. Therefore, the limit in (35) is 
well defined and provides us with a finite answer for the quantum correction to the kink mass. Despite 
appearances, formula (35) gives the same results as the CCG and DHN formulas, as one can check from 
standard Tables of special functions. 



3.1.4 The A'-heat kernel asymptotics 

The main lesson of our previous calculations is that the a = iV-kinks are very special because the kink 
wells are transparent. This fortunate fact allows the computation of the kink semiclassical mass by 
means of relatively simple formulas extracted from the fairly complete spectral information about the 
corresponding Schrodinger operators. Dealing with more exotic kinks and self-dual vortices of various 
types, see [27] and [18] for reviews, a method has been developed that only requires the knowledge 
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of the potential well V{x) in the computation of one-loop mass shifts of low dimensional topological 
defects. Starting from the Gilkey-DeWitt-Avramidi high-temperature expansion of the heat function 
kernel [23, 24, 25, 28, 29], the idea is to express, via the Mellin transform, the energy shifts AE(<Pk) 
as a truncated series whose coefficients are the Seeley coefficients of the K-heat function. We stress 
again that this strategy only requires knowledge of the kink potential well V(x), which encodes all the 
spectral information; the Seeley coefficients are integrals over the real line of the Seeley densities, built 
from polynomials in V(x) and its derivatives. Therefore, the virtue of this procedure is the universal 
application to any (l+l)-dimensional one-component scalar field theory model supporting kinks, see [19] 
and [27]. The asymptotic approach, however, requires arduous calculations which can be alleviated by 
using a symbolic software platform, although it has been successfully generalized to (l+l)-dimensional 
iV-component scalar field theory [8, 30, 18] and even to (l+2)-dimensional Abelian Higgs models [31, 32]. 
Briefly, the key formula giving the one-loop kink mass shift has the form of a truncated series 3 : 

1 Nt 1 
8vr ^ 20^0 

where 7(0, b] is the incomplete gamma function, c n (K) are the associated Seeley coefficients of the K-heat 
function expansion, Nt € N is the truncation order of the series, and (3q is a parameter that freely sets 
the integration range in the Mellin transform in such a way that the outcome will be optimized. We test 
the method by applying this formula to the a = iV-kinks, a case where we know the exact result. In fact, 
because we know the K^-heat and Kq exactly we merely expand (34) in powers of (3 

-N 2 (3 00 nn+1 ( N \ . p -(3v 2 00 

'"' ""' . < 37) 

using the very well known series expansion of the Error function, to immediately identify the Seeley 
coefficients : 

Plugging these coefficients (38) into expression (36) and taking large enough Nt and /3o, we reproduce 
the results for the one-loop kink mass shifts given in Table 2. 

In general there is no exact information on the if-heat function. Nevertheless, the Seeley coefficients 
c n {K) can be identified from the if- heat kernel expansion: 

c n (K)= / dxc n (x,x) = / dx {0) C n {x) . (39) 



The Seeley densities c n (x,x) can be obtained by solving the recurrence relations -coming from the series 
expansion of the fundamental solution of the if-heat equation- 



L '/A iPV 



Wc -w-^[ Mc -w-gG) 

with "initial" (infinite temperature) conditions: ^Cq(x) = 5 k0 . 



dxi 



^ r, 



(40) 



3 A similar formula also works for planar self-dual vortices. 

4 Note that in formula (36) we do not need the first two coefficient: co(K) and c\{K). We have shown in previous 
works that the contribution of Co is removed by zero-point renormalization whereas the ci term disappears after mass 
renormalization. 
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In order to test the asymptotic approach we now compute (36) from the coefficients obtained as 
solutions of the recurrence relations (39) and (40). We shall choose a truncation order Nt = 40, e.g., in 
the fourth member TV = 4 of the hierarchy: U^{(f>) = |[1 — 2cos(| arcsin 0)] 4 , \<f>\ < 1. The Seeley 
coefficients arising from the recurrence relations are shown in Table 3. 



n 






n 






n 






n 




1 


40.0000 




11 


1313320.5922 




21 


1546743.9642 




31 


12813.7015 


2 


266.6667 




12 


1825331.1803 




22 


1151061.4774 




32 


6508.5467 


3 


1386.6667 




13 


2335057.9936 




23 


818531.0868 




33 


3204.2076 


4 


5699.0476 




14 


2766565.6064 




24 


557297.1794 




34 


1530.3678 


5 


19121.4392 




15 


3052196.9105 




25 


363948.9646 




35 


709.7358 


6 


53853.5835 




16 


3150326.7363 




26 


228360.0592 




36 


319.8809 


7 


130166.9684 




17 


3054683.3061 




27 


137877.7460 




37 


140.2218 


8 


274845.7626 




18 


2792761.2600 




28 


80219.7711 




38 


59.8280 


9 


514360.2887 




19 


2415316.3095 




29 


45035.6583 




39 


24.8636 


10 


863458.4292 




20 


1981777.3296 




30 


24426.1190 




40 


10.0713 



Table 3: Seeley coefficients of the K^-heat function expansion. 



These data coincide exactly with those obtained from the analytical formula (38) for N = 4: 

nn+l 

Cn(K(4)) = (2 n-l)H (l2 "^ + 22 " _1 + + 42 " _1) • 

The 40 Seeley coefficients, obtained either way, are now plugged into formula (36). The tuning of the (3q 
parameter to the value /3q = 1.3125 offers an optimum estimation of the one-loop kink mass shift 

AE(4$) -1.58003 

to be compared, e.g., with the answer given by the CCG formula. The link between N-t and (3$ can be 
understood by looking at Figures 3(b) and 3(c): there, several partial sums, as well as the total series, are 
plotted as functions of /3. For each iV^-order there is a point (3q in the /3-half-line where the partial sum 
is closer to the exact value. The further we can push f3o towards infinity, the larger part of the function 
we can pick in the integral, and the optimum value is reached. 

3.2 The a £ N case 

In this section we shall apply the generalized DHN formula and the asymptotic approach to compute the 
one-loop (15)-kink mass shifts when a ^ N. In this case A a \k), given in (23), are not null and, moreover, 
there is no half-bound state in the spectrum of K^ a \ 

3.2.1 The generalized DHN formula 

In [8] the extension of formula (31) to this more general situation has been described thoroughly. The 
"generalized" DHN formula reads: 

AB( te ) = ig^ + i w -H + ±jr*^>^i^ + i(v(.))-i(vw) <«> 

where w| are the bound-state eigenvalues of K = — ^ + v 2 + V(x) and si is equal to | if the higher 
eigenvalue ojf = v 2 lies in the continuous spectrum threshold, but it is otherwise one: u? < v 2 . The new 
ingredient in the formula encodes the imbalance between the half-bound states of K and Kq. 
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The application of formula (41) to the a ^ N models only requires knowledge of the spectral density 
derived from (23) and the consequent total phase shift: 

dS [a] .,x (/ n (/ , N ij(ik)T 2 [ik}T 2 [l - ik] - ^(1 - ik)T 2 [-a}T 2 [l + a] 

il){-ik)-i)(-<J-ik)-^{l + a-ik)+ Vy ' i, J -M^r-n Ur tWf, V 1 • ( 42 ) 



r/A- 



r 2 [i - ik]r 2 [ik] - r 2 [-a]r 2 [i + a] 



Plugging (42) into (41) for the a N-kinks we obtain: 

a{cr+ 1 



1 I[o] 



3 - 



j=0 

Owing to the complexity of the 
possible. Thus, we extract the 
results: 



2vr 



a 1 

4 + 27 



dk 



'o 



dk 



Vk 2 + a 2 + ±= ; 



<7" 



derivative of the phase shifts, exact integration of this formula is not 
information by means of numerical integration, offering the following 
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0.0 


0.00000 




1.0 
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One-loop parent kink mass shifts for several values of a £ [0.0, 4.0] estimated by applyin 



generalized DHN formula. 




Figure 4: Graphical representation of: (a) one- loop kink mass shifts in the range: a G [0.0,4.0]. (b) Seeley 
coefficients within the same range. 



The values depicted in Table 4 have been drawn in Figure 4 as solid points. The figure shows a 
remarkable pattern: the continuity of A-E^^] as a function of the parameter a. The generalized DHN 
formula (derived in the mode-number cut-off regularization procedure plus a ^-weight on the half-bound 
states coming from the ID Levinson theorem) guarantees this behaviour. Had we applied the original 
DHN formula for all values of a € M we would have found discontinuities at integers values of a. Note 
also that AE[4>^]/mdh tends to zero when a vanishes, as should be the case, see Figure 4(a). 
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3.2.2 The asymptotic approach 



Finally, we shall apply the asymptotic approach in the a £ N models. In this case, computation of the 
Seeley coefficients c n (K^) in (36) can only be achieved by solving the recurrence relations (40). In 
Reference [20] a symbolic program code has been developed that performs this laborious task. We use 
this code to calculate the Seeley coefficients, adapted to the present models. A graphical representation 
of these coefficients is shown in Figure 4(b) for different values of the a parameter with a separation of 0.1 
in the range [0.4, 2.0]. In general, the larger the value of a, the larger the coefficients c n (K^). Plugging 
these coefficients into formula (36) we estimate the quantum correction. Frequently, the response of (36) 
has been evaluated by employing a series involving ten terms, i.e., Nt = 10, although this can vary along 
the values of the parameter a. The results obtained are shown in Table 5. 
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Table 5: One-loop (15) kink mass shifts for several values of a G [0.0,4.0] estimated by applying the 
asymptotic approach. 

The values seen in the Table 5 are drawn in Figure 4(a) as empty squares in order to appreciate the 
concordance between the data obtained upon using the generalized DHN formula and the asymptotic 
approach. Only a small difference between the data arises for a < 0.6 and lower. For a small the series 
in (36) converges slowly because of the factor (cr 2 ) 1- ™ in each summand; a larger number of terms is 
needed to obtain a precise answer. The limitations in the computations of the Seeley coefficients make 
the asymptotic approach less exact in the range of values of the parameter a < 0.4. 

4 Conclusions and further comments 

We have constructed a one-parametric family of (l+l)-dimensional one-component scalar field theory 
models governed by the action (2) with the potential (18). This family includes two famous members: 
the sine-Gordon and 4 models. The models in this family are characterized by very specific kink wells: 
the second-order kink fluctuation operators are Schrodinger operators with Poschl- Teller potential wells. 
In general, the family of potentials (18) is only of C 2 (M) class as a function of the field </>: Derivatives of 
order higher than two are discontinuous or divergent at the vacuum points. There exists a similar family 
of models, with class C°° potentials, for two or more scalar fields. For example, the (l+l)-dimensional 
two-component scalar field theory model with a C°°(]R 2 ) potential 

U{<PiA2) = + 2a4> 2 2 - l) 2 + 2a 2 4>j4 
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admits the kink solution 

(j)f(x) = - tanhx , cpf(x) = 

whose second-order small fluctuation operator is: 



K {*) = ( -i& + 4-6sech 2 x ^ ^ 

dx 2 



—f? + a 2 -a(a + l)sech 2 z 



Here, the generic Poschl- Teller potential well arises in the study of the orthogonal fluctuations to the 
kink solution, see [8]. 

Finally, the spectral problems associated with kink stability analysis in the models studied in this 
paper are solvable and therefore the one-loop quantum correction to the kink mass in these cases can 
be computed exactly. By the same token as before, AE^-[(j)j? '], the one-loop kink mass shift due to 
orthogonal fluctuations to the kink, is a continuous function of a. Note, however, that our methods have 
been improved such that the agreement between the results via the generalized DHN formula and the 
asymptotic method is now much better than eight years ago, see [8]. 
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